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1/^ ' Abstract 

In this paper, the symmetry group of a differential system of n quadratic 
r ^ ] homogeneous first order ODEs of n variables is studied. For this purpose, 

we consider the action of both point and contact transformations to signify 
the corresponding Lie algebras. We also find the independent differential 
invariants of these actions. 
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^ '■ 1 Introduction 

QQ , A quadratic homogeneous differential system (QHDS) is a generalization of sec- 

f^ ' end degree polynomial system of differential equations that is rooted in natural 

r — I , science, and also in mathematics by virtue of Hilbert's sixteen problem of his 

^^ ■ famous list of problems (See, e.g., [T] and references therein). They are studied 

00 , in some different view point of their applications. For example, in physical sci- 

^^ ' ences, in van der Pol oscillator as an important example of qualitative theory of 

ordinary differential equations (ODEs), in mathematical echology and in partic- 
ular in Volterra-Lotka equations, and many other applications in astrophysics 
and fluid mechanics [l]. 
C^ ' There are many aspects for studying a QHDS. For instance, in [3] geomet- 

ric classification of trajectories of a given two dimensional QHDS has studied 
for determining the invariant lines through the origin of the QHDS and also 
its location in the plane of parameters for classification of the geometry of its 
trajectories. In [8], the necessary conditions for the existence of polynomial first 
integrals and a polynomial symmetry field of a n-dimensional QHDS have given. 
At the present work, the symmetry group of a QHDS of every finite di- 
mension is studied. We specify symmetry groups of both point and contact 
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transformation groups. 

Assume that we have a differential system £ oi n quadratic homogeneous 
first order ODEs of n variables, in which each of them is a function of t: 

£ : ii^ Fi{xi,- ■■ ,Xn), Fi = '^a)kXjXk, (1.1) 

where the coefficients, a^, ^ for 1 < i,j, fc < n, are arbitrary real constants. 
The solution space of differential system (II. ip is an immersed submanifold of 
the first order jet space, J-'^(R, R") of dimension 2n + 1 and local coordinates 
(t,x, p) := (i, a;i,-- ■ ,a;n,pi, • • • ,Pn), when each pj is just Xj. 

According to [6] , the action of contact transformations group is introduced 
on the jet space J^(R, R"), and any element of this group is as the following 
transformation: 

r = x(t,x,p), X, = i/;,(t,x,p), Pj=^,(t,x,p), (1.2) 

where 1 < «, j < n. An important special case of contact transformation group 
is point transformation group, where the action reduces to the following change 
of coordinates 

r = x(t,x), X, = v^,(x,p), 

when the variables come from the jet space J'^(R,M"). 

Remark 1.1 In general, the point and contact symmetry groups of a geomet- 
ric object are not necessarily equal. A jet space J'^(R™,R") of order k can be 
imbedded in other jet spaces j'(R'",R") of order I > k as a submanifold. Also 
a contact transformation o/ j'(R™,R") can be restricted to j'^(R'",R") and 
then acts as a point transformation. Thus, the point transformation group of 
an geometric object is a Lie group of the contact transformation group and in 
the special case of symmetries, we conclude that the point symmetry group is a 
Lie subgroup of the contact symmetry group. 

In the next section, we will see that point symmetry group of the differen- 
tial system £ is of dimension 1, where its contact symmetry group has infinite 
dimension. 



2 Point Symmetry of a QHDS 

As is indicated above, an important question of a quadratic homogeneous sys- 
tem is about its symmetry group. There are various methods for signifying 
symmetry group based on Lie's method of infinitesimal generators, that some 
of them are exist in [2] , [4] , [5] and [6] . In this work, we will follow the method 
of [6] for finding point and contact symmetries. 



In jet space J^(R, K") the system of equations (|l.ip will be as follows 

Pi = ^ a)k Xj Xk, ^ <i,j,k < n. (2.3) 

For finding symmetry group of system £, we consider the general form of the 
infinitesimal generators provided from a point transformation as follows: 

i—l 

in which, T and KiS are smooth functions of variables t,xi, ■ ■ ■ ,Xn. 

The fist order prolongation of v will be then as the following expression: 



d_ 
dpi' 



^^'^-- + E^*£' (2-5) 



where for each 1 < i < n we have 

Kl = Ku, + J2 K^^o., - P^{Tt + Y. ^^-.)- (2-6) 

J 3 

According to ^, v^^' as a prolongation of w is a symmetry of system £, if and 
only if it satisfies in relation w'^^^[f] — 0. By effecting v'^'^^ on the equations of 
the system we conclude the following system of equations for i = 1, • • • , n: 

Kl ( Y. («U- + O^klW) +---+Kn[Y. «fc + «fcn)^fc) 
k k 

-Ki^t - Pi Ki^x^ - p2 Ki,^^ pn Ki^^^ (2.7) 

+Pi{Tt + Pi r^Ei H h Pn Tx„) = 0. 

In the last equation, since t,Xi,pj (1 < i,j < n) are arbitrary and functions Ki 
and Nj only depend on i,Xi, so equations (|2.7p will be satisfied if and only if 
we have the following system of equations 

k k 

X,,,, = 0, K,^^, = 0, • ■ • , X,,,,_, = 0, X,,,.^, = 0, • ■ • , X,,,„ = 0, (2.9) 

Tt-if.,,, =0, T,, =0, ■•■, r,„=0. (2.10) 

for every 1 < i < n. At first, from Eqs. (|2.9p - (|2.10p we understand that T 
depends only to t and for each i, Ki is a function of t and Xi. For a fixed 
1 < i < n, by replacing the Eq. Tt = Ki^^i in (|2.8p and then differentiating with 
respect to 0:^, we will find the following expression 



rt(^(a^fe+aL)a;fe)=r„, 



By solving the last equation in respect to t, we have the following solution 
Tit) = ci exp (t Y^ [a],, + <,)xfc) + ca, 

k 

for constants ci and C2. Since T is just a function of t, so ci must be zero, and 
therefore 

1 [Lj Xi, ■ ■ ■ , Xn) — C, 

for arbitrary constant c. Then for each i, we see that i^Tij^. = and so is 
dependent to t only. But by solving Eqs. (|2.8p with respect to Ki{t) we obtain 
the below solutions when i varies between 1 and n: 

K,{t) = exp(i^(a^fc+4Ja;fe) 

k 

ljexp(^-tJ2 {<k + <d^k) Fdt + cA, (2.11) 

where we assumed that F = J2j=^i ^j (^) ( J2k ('^jfc + ^kj )^k ) and C3 be arbitrary 
constant. Eqs. (|2.1ip and the fact that all of KiS are just depend to t implies 
that 

Ki{t, xi, • • ■ , Xn) — for 1 < i < n. 

when c is an arbitrary constant. The symmetry group of the general infinitesimal 
generators, that provided as w = c-^r will be just the translation of t by a constant 
coefhcient of the parameter s of the 1-parameter subgroup: 

{t,Xi,-- ■,Xn) ^ icS + t,Xi,-- ■ ,Xn) 

Theorem 2.1 Set of all infinitesimal generators of the point symmetry group 
of the system il.l]) is a 1-dimensional Lie algebra with the base {^}- There- 
fore, the symmetry group of the system \1.1]) is a 1-dimensional Lie group of 
time translations. 

According to theorem 2.74 of [6], the invariants u = L{t,xi,X2, ■ ■ ■ ,Xn) of 
one-parameter group with infinitesimal generators of the form (|2.4p satisfy the 
linear, homogeneous partial differential equations of first order: 

v[L] = 0. 

The solutions of the last equation, are found by the method of characteristics 
(See [6] and 4 for details). So we can replace the last equation by the following 
characteristic system of ordinary differential equations 

dt dxi dx2 dxn 



By solving the equations (|2.12p of the differential generator v — c-^ , we (locally) 
find the following general solutions 

Il{t,Xi,X2,-- ■,Xn) = Xi = di, •••, In{t,Xi,X2,--- ,Xn) =Xn = dn, (2.13) 

for constants di when i — 1, • • • , n. 



Theorem 2.2 The independent first integrals of the characteristic system of 

.A 



the infinitesimal generator v = c4^ are the derived invariants V2.13\) . 



3 Contact Symmetry of a QHDS 

On the other hand, we consider that the infinitesimal generators comes from 
contact transformations, that is, the contact transformation group is acting 
on the jet space J^. Hence, we may suppose general form of a infinitesimal 
generator as 

-^Ih-Sa-.A.Sp.A. (,.14) 

i=l * i = l ^■' 

for arbitrary smooth functions T, Ki, Pj of variables t, Xi,pj. 

Since this expression depends to variables of 1-jet space and so is a infinites- 
imal generator of the jet space, so it dose not need to be prolonged. Therefore, 
according to [3] f is a symmetry oi £, if and only if, it satisfies in v[£] ~ 0. 
After acting v on the system we obtain the relation 

Ki Y^ (a'ifc + a.l^)xk + --- + KnJ2 i<k + a^kJxk - P^ ^ 

k k 

where i varies between 1 and n. Thus the general form of each symmetry of the 
system will be as follows 

1=1 j,k 

in which T, Ki, K^, ■ ■ ■ , Kn are arbitrary functions of t, xi, ^2, ■ • ■ , a;„. One may 
divides v to the following vector fields 






d ■s-^ , , , , d 






Vn+1 = K^{^^+J2iai, + aiJxk^), 

j,k J 



(3.16) 



Lie bracket of each two of wi , • • • , Vn+i is defined to be the commutator [vi ,Vj] := 
Vi Vj — Vj Vi {1 < i, j < n + 1). It is not hard to see that the commutator of 
any two of t;i, • ■ • , f„+i is a linear combination of them, and hence vi, ■ ■ ■ , Vn+i 
generate a Lie algebra. For instance for 1 < a, (3 < n ii we assume that 



dXa 



j,k 



vp+1 



K. 



d 






then we have 



d 

j,k 
j k 

k 



T-klil^k 



kf3 



IXk 



-^P ( Yl ^"^P' ("m/3 + 4m)^") ( Yl (°«fe + «ia)2^fe) } Q- ■ 



l,m 



One can compute the table of Lie symmetry algebra of commutators arising 
from infinitesimal operators p.l6p that is given in table 1. 



Vl 



V2 



Vl 
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V2 


~Vl 


V2 


Vl 


-V2 







Vn 


Vn 


~Vi 


V2 


-Vn 


Vn+1 


Vl - 


w«+i 


V2 - 


- W„+l 



Vn 



Vn+1 



Vn - Vl W„+i - Vl 

Vn - V2 W«+l - V2 

U„+i - Vn 
Vn - Vn+1 



Table 1. The commutators table of Lie algebra of equation (|3.14[) 
Therefore we have the following theorem: 

Theorem 3.1 The symmetry group of the system £ of quadratic homogeneous 
differential equations, \1.1\) . under the action of contact transformations is an 



infinite dimensional Lie algebra generated by infinitesimal operators 113.16]) . 

By referring again to theorem 2.74 of [S], and considering the infinitesimal 
operators of the form (|3.15p . the invariants /(t, x, p) of their one-parameter 
group, are signified by the characteristic system of ordinary differential equations 
as follows (i, J, /c = 1, • ■ • , n) 

M ^doH ^ dpj 

T K, K, Z, (a{, + ai,)x, + ■ ■ ■ + K^Y.k Kk + O^ 

in which determines 2n invariant functions. By solving these equations we find 
that 

/i(i,x,p) = I {Kidt-Tdki)^di, 

/„(i,x,p) ^ f{Kndt-Tdkn)=dn, (3.18) 

In+i{t, X, p) = / ( X! ^J ("ife + ^h )^fe) ^* ~ Tdpi^ dn+1 



hn (i, X, p) = / ( X! ^^ ("ife + "fcJ )xkjdt- IT dpn = d2n 

Theorem 3.2 The derived invariants Ii3.18\) as independent first integrals of 
the characteristic system of the infinitesimal generator Ii3.15\) . provide the gen- 
eral solution 

S{t, X, p) := ifiih {t, X, p), hit, X, p), • • • , l2„it, X, p)), 

with an arbitrary function Lp, which satisfies in the equation v[ip\ — and there- 
fore in 

1=1 j,k 

The solution of the last equation is defined implicitly by S{t, x, p) = 0. If ^ ^ 0, 
then the solution can be written explicitly by t — 7(x, p). 



Example 3.3 Clearly, by assuming T = 1 and Ki = we return back to a 
point symmetry, in which as we saw in previous section, the symmetry group 
provides by translating time and fixing other variables. 



Example 3.4 If we suppose T — 1 and Ki = 1 for some i and Kj = for 

j =f^ i, then we have the foUowing expression of infinitesimal generators 

Its one-parameter group then will be as follows 

(i,x,p) H^ (^t + s,xi,---,x^ + s,---,Xn,Pm + a'^s^ + s^{a'l'^ + a]^i)xkj 

k 

where s is the parameter of the flow and m changes over 1, • ■ • , n. If the point 
(i, X, p) be fixed, then the fiow of w, by restricting the jet space coordinates to 
the i, Xi and pj axis, for a 1 < j < n and the fixed value i, is a space parabolic. 

Example 3.5 Let T = t, Ki = 1 for some index i and Kj = for other 
indices j ^ i. The infinitesimal operator ti and its flow are given respectively 
by changing the coefficient of ^ to i in (|3.20[) and the component t + s of right 
hand side of dS^Ol) toie". 

By fixing (t,x, p) and acting the symmetry group on it, then the flow of u, 
under the projection to the coordinates in terms of i, x,; and Pj for a 1 < j < n, 
will be similar to the space curve (s, e'', s^ + s) and has two branches, the first 
of space exponential map and the second of a space parabolic. 

Example 3.6 For the case which we assume that T = ct when c is a constant, 
and Ki = 1 when i = 1, ■ • ■ , n, the contact infinitesimal operator v is as follows 

I j,k 

The one-parameter of v then is 

(t,x,p) ^ (te^^x + s,p,„ + s2^a;7 + s^(a™+aS)a;fe) 

I l,k 

when in it, m varies from 1 to n and we assumed that 1 < fc, / < n. The action 
of the symmetry group on a fixed point, will give a flow that its projection to 
three arbitrary axis {t,Xi,pj) is similar to the 3-dimensional flow explained in 
example 3.5. In the particular case of this example, the characteristic equations 
tend to the following independent differential invariants for 1 < a < n and 
n+l< f3 <2n: 

/a(i,x,p) ==i(l -c/ca), If3{t,yi,p) ^ t i^J2j,kia^jk'' + akj'')xkj -ctpp^n- 



Example 3.7 By assuming T — and Ki = Xi we find that the infinitesimal 
operator as 

I j,k 

and the one-parameter group of v as 

(t,x,p) H^ (t,e"x,p, + (e2''-l)^(a}fc+<j-)xjXfc 

where i varies on values 1, • • • , n. Its flow also signifies a parabolic if we reduce 
the jet coordinate (t, x, p) to (t,xi,pk) for some specified I and k. The concluded 
invariants then will be as follows (1 < ct, /3 < n) 

/i = t, la = IuXq, — lna;Q_i, 

Finally, we present an example of the case in which KiS arc functions dependent 
to pjS. resp. 

Example 3.8 As wc indicate, we suppose that T ~ and Ki = pi for i = 
1, • • • , n. So the infinitesimal generator is in the form of 



^-ll{p^^ + [llpM)k+ai,)x,)^V 



dpi 

that has the following trajectory in term of the parameter s 

(i,x,p) ^ (t,x + ps,p') 
where p' consists of n component, which its z*'' component is 



Pt^Pi+ [ cxp 



(E(«^fc+«HK)5)-l)[- 



T,Mk + (^lt)^k 

One can derive the following independent differential invariants 

il — t, ia — PaXa—l Pa — lXaj 



for a, (3 — 1, • • • , n. Clearly, projection of the jet coordinate to the coordinate 
(xmiPi) for specified amounts m and I, is the graph of the exponential map. 



According to [7], if G acts smoothly and transitively on a manifold M, then 
M is isomorphic to G/H as a homogeneous space, which obtained by quotienting 
by a closed Lie subgroup that is, an isotropy subgroup provided by element of 
G. Thus, we can use of the exponential map, exp : TM — > Af , defined in a small 
neighborhood of identity element of M as a Lie group. From Baker-Campbell- 
Hausdorff formula we know that for vector fields X, Y of the Lie algebra of a 
Lie group, if we have [X, Y] = 0, then the flow oi X + Y \s £,t = <t>t°'4't = 4't°'Pt 
where (f>t and ipt are flows of resp. X and Y. 

Remark 3.8 With similar method as indicated in above examples, by differ- 
ent selections of coefficients of i3.15\) as the linear combinations of 1-jet space 
variables, their corresponding invariants are in the forms of i3.18\) . and their 
proper projections of trajectories to 2 or 3-dimensional space of their variables, 
providing their commutator be zero, are similar to parabolic, exponential map 
or compositions of two branches in the similar shape to them. 

This result is inferred, since the symmetric group prepared in the above 
examples act transitively, if we add the condition that the Lie bracket of every 
two infinitesimal operators be zero, then the fiow of every linear combination of 
any two of infinitesimal generators is a composition of fiows made by these two 
operators. 
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